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dynamics with five or six active flavours. The new results can be used to study 
the stability of a s evaluated at a high scale from the knowledge of its value at Mz- 
We furthermore derive a low-energy theorem which allows the calculation of the 
coefficient function of the effective Higgs boson-gluon operator from the decoupling 
constant. This constitutes the first independent check of the matching coefficient 
to three loops. 
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1 Introduction 



The decoupling of particles with masses much heavier than the considered energy scale 
has a long history [TJ. It is tightly connected to the construction of an effective theory 
containing only the light active degrees of freedom in the dynamical part of the Lagrange 
density. Within the framework of QCD decoupling constants for the strong coupling 
a s are known at two- [2HI], three [I] and even four-loop order [5j[6]. Recently also the 
expression for the simultaneous decoupling of two heavy quarks has been computed at 
the three-loop level [TJ. 

Decoupling relations are also important in the context of supersymmetry where the Stan- 
dard Model constitutes the effective theory. Two-loop corrections for a degenerate super- 
symmetric mass spectrum are known from Ref. [8],[9] and the general result can be found 
in Ref. [TO]. In this paper we compute the three- loop corrections for several different 
assumptions on the masses of the MSSM. 

There is an interesting connection between the decoupling constants and the effective 
coupling of a CP neutral Higgs boson to gluons which is defined via the Lagrange density 
(the superscript marks bare quantities) 

£y, c ff = ~^C%^ + £qcd> W 

with 

0\ = \g%G ^. (2) 

where <fi is the Higgs field v is the vacuum expectation value and the field strength 
tensor in QCD. Cq CD is the QCD Lagrange density with five active flavours. The first 
term in Eq. (QQ) describes the coupling of the Higgs boson to two, three and four gluons. 

In Ref. [I] a low-energy theorem (LET) has been derived which connects C\ to the deriva- 
tive of the decoupling constant for a s with respect to the top quark mass. As far as 
supersymmetry is concerned a next-to-leading order (NLO) version of the LET has been 
derived in Ref. [11] (see also Ref. [12]). In this way the NLO supersymmetric QCD 
(SQCD) corrections to C\ obtained in Ref. [13] could be confirmed. We re-derive the 
LET, apply it at three loops and thus obtain the coefficient function C\ which is needed 
for NNLO prediction of Higgs boson production and decay within the MSSM. With our 
calculation we confirm the result for C\ obtained in Ref. [T^fTo] by an explicit calculation 
of the vertex diagrams. 

The outline of this paper is as follows: In the next Section we describe the calculation of 
the decoupling constant for a s to three loops and discuss the numerical influence in the 
computation of « s (Mgut)- Afterwards we derive in Section [3] an all-order low-energy- 
theorem which we use to compute C\ to NNLO accuracy. We summarize and conclude 
the paper in Section HJ In the Appendix we present a compact expression of the exact 
two-loop result for the decoupling coefficient. 
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2 Decoupling of heavy super symmetric particles 



In order to compute the decoupling effects of heavy particles from the running of a s 
one can use the well-established formalism derived in Ref. jl]. It has been applied to 
supersymmetry in Refs. [HHTU] where two-loop corrections have been computed. 

The starting point is the relation between the strong coupling in the full theory, which is 
in our case the MSSM, respectively, SQCD, and the effective theory, QCD 

«i QCD) (/i) = C as (/i)«f QCD) (/^)- (3) 
At that point some comments are in order: 

• ai^ 013 ^) is defined in the five or six flavour theory, depending on whether the top 
quark is integrated out together with the supersymmetric particles or not. 

• a^ C (jtz) is defined in the MS scheme based on Dimensional Regularization 
(DREG). a^ CD \fi) is defined in the DR scheme since the supersymmetric the- 
ory is regularized using Dimensional Reduction (DRED). DRED is implemented 
with e scalars, where the details can be found in Refs. [T5l[T6l. 



Ca s (/-0 as introduced in Eq. (j3]) has two tasks: (i) it has to decouple the heavy 
particles not present in the effective theory, and (ii) ( ag (/i) has to ensure the change 
of regularization from DRED to DREG. In principle the two tasks can be performed 
in two steps as it has been proposed in Refs. [5HTU] . However, it is more convenient to 
choose the same renormalization scale for the decoupling and the change of scheme. 
Calculations along these lines have also been performed in Ref. [T^[T7]. 

In principle each vertex containing a s can be used in order to compute ( ag . It is, 
however, convenient to use the gluon-ghost vertex in order to compute the decou- 
pling constant via jl] 



where the superscript "0" marks bare quantities. 0, £3 and (3 are the decoupling 
constants of the gluon-ghost vertex, ghost and gluon propagator, respectively. They 
are obtained from the hard part of the corresponding Green's function. The corre- 
sponding formulae can be found in Ref. jl] where a derivation has been performed 
in the framework of QCD. It can be taken over to SQCD without modifications. 
The renormalized decoupling constant is obtained from 

s 

where Z as and Z a i a are the renormalization constants for a s in the full and effective 
theory, respectively. 
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All occurring parameters are renormalized in the DR scheme, except the e scalar 
mass which is renormalized on-shell with the condition M £ = 0. The corresponding 
counterterms can, e.g., be found in Ref. [TJ~ 



Assuming a strong hierarchy among the quarks one encounters in the case of QCD vac- 
uum diagrams which contain only one mass scale. The occurring integrals can even be 
computed up to four-loop order [HIE]. Two scales appear if two quarks are integrated out 
simultaneously. This has been done in Ref. [7] to three-loop accuracy. 

In the case of supersymmetry significantly more mass scales have to be considered. In our 
approach we have the gluino and top squark masses (rrig, m^, m t - 2 ) and a generic squark 
mass rrtq which we take as the average of the up, down, strange, charm and bottom 
squarks. In addition there is the e scalar (M e ) and the top quark (m t ) mass. The latter 
only appears if we match to five-flavour QCD since m t = is chosen for the matching to 
six-flavour QCD. Up to two loops ( ag can nevertheless be computed exactly [10] taking 
into account the dependence on all mass parameters. The analytical result can be found 
in the Appendix. At three-loop order, however, approximations have to be adopted in 
order to be able to compute the integrals. Motivated by scenarios which are currently 
discussed in the literature we have chosen 
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(h2) 
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> m t , 


(h3) 


nig p 
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ps m g 


> m {l 


~ m t , 



(6) 

where in the case of "^>" an asymptotic expansion in the corresponding hierarchy is 
performed. In the case of "~" a naive Taylor expansion in the difference of the particle 
masses is sufficient. For all hierarchies we assume that M £ is not zero but much smaller 
than all other masses. In this way we ensure that the e scalar is integrated out and not 
present in the effective theory. Thus, in the latter dimensional regularization can be used. 
In what follows the heavy mass scales for each hierarchy are also denoted by msusY in 
case they are identified. 

At three-loop order terms up to 0(l/m^j SY ) have been computed for (hi) and (h3) and 
up to G(l/m~ ) and 0(l/m^ VSY ) for (h2). For each mass difference at least four expansion 
terms (i.e. terms including {mf — m'j) 3 ) could be evaluated. It is either possible to expand 
in the linear or the quadratic mass difference. Formally both choices are equivalent, 
however, in practice it turns out that depending on the actual numerical values of the 
parameters one can be significantly better behaved than the other. Similarly there is 
a freedom to choose a mass parameter, m R , around which the expansion is performed. 
m R should be of the order of the involved masses. Note that for (hi) and (h2) only one 
reference mass ijir is required whereas for (h3) one needs two as can be seen from Eq. ()6]). 
Again there may be significant numerical differences and thus we adopt the following 
choices when evaluating the three-loop corrections to the decoupling coefficient 



m fl + mi 2 + 10m g + nig 
- — 



(hi) m R = m- t , m R = m~ t ,m R = m g ,m R = m q ,m R 
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(h2) m R = m i2 ,m R = m g ,m R = nig ,m R 

(h3) m Rl = m- t2 , m Rl = m g , m Rl = rrig , m Rl ■ 

m t + m~ t 

m R2 = m t , m R2 = m~ tl , m R2 = 



m~ t2 + lOrrtg + m-g 

12 ' 
mi 2 + 10mg + m g 

12 



(7) 



In the following it is convenient to consider the perturbative expansion of ( ag which we 
define as 
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where the /x dependence of ai SQCD ' ) and Q'i is suppressed on the right-hand side. 

The general results are quite lengthy and will not be presented in this paper. However, 
in order to get an impression of the results we present ( as for the hierarchy (hi) with a 
degenerate supersymmetric mass spectrum which reads 
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(9) 



where x t s = m t/ m svsYi h — hi(/i 2 /mf ), Is = ln(/i 2 /mg USY ) and l x = \n(x t s)- The ellipses 
denote terms of order xfg. The corresponding results where the matching is performed 
to six-flavour QCD, i.e. where the top quark is not integrated out and thus treated as 
massless in the loop integrals, reads 



C (1) = ---la 
>(2) = ZZ _ _Z_ Z . Z 2 



s ■ 



5 



, / o\ 11203 1495, 541 ,„ 4657 A/ , , , 

® = - 4608 - ™r' s+ 288^-^ + 92l6 C(3) - (10) 

All analytical expressions corresponding to the hierarchies of Eq. (jHJ) can be found in the 
file decsusy31.m obtained from [19]. 

Let us in the following test our approximation at two loops by comparing to the exact 
result. For this purpose we adopt the following values for the input parameters 

m t = 150 GeV , A t = 100 GeV , Mq s = 500 GeV , /i SU sy = 100 GeV , 
tan/3 = 10, M z — 91.2 GeV, sin 2 6 W = 0.2233 . (11) 

where A t is the trilinear coupling, /^susy is the Higgs-Higgsino bilinear coupling from the 
super potential, tan/3 is the ratio of the vacuum expectation values of the two Higgs 
doublets, Mz is the Z boson mass, 0w the weak mixing angle and Mq 3 , a soft SUSY 
breaking parameter for the squark doublet of the third family. Furthermore we set the 
renormalization scale to fi = 500 GeV. These parameters can be used to compute m^, 
m~ t2 and 9 t as a function of the singlet soft SUSY breaking parameter of the top squark, 
Mu 3i r (see, e.g., Ref. [20]) by diagonalizing the corresponding mass matrix. The result is 
shown in Fig. Ufa). Furthermore we choose for simplicity mj = rrig = trig. This allows us 
to consider in Fig. H^b) both the exact result for Q a (solid line) and the approximations 
(dahed lines) based on the hierarchies (hi), (h2) and (h3). The latter are obtained from 
the (naive) averages over the various representations, i.e., the different choices of 
according to Eq. (JTj). One observes that in the whole range of Mu 3j r at least one of the 
hierarchies approximates the exact to a high degree, which provides the motivation to 
proceed in a similar way at three loops. 

Since at three-loop order the exact result is not known a criterion is needed in order to 
select the best approximation among the various choices at hand. For this reason we 
define 



<5app 



A2) _A2) 
Sapp Sexact 



^ exact 



A3)c _ A3) 
Sapp S>app 



A3) 
Sapp 



(12) 



where "app" marks an approximation result and the superscript "c" indicates that the 
highest terms in the expansions are cut. For each set of input parameters we choose the 
representation which leads to the minimal value of <5 app . The first term on the right-hand 
side of Eq. (|12|) guarantees that the approximation works well at two-loop order whereas 
the second term assures the convergence of the expansion. 

The three-loop result (a] is shown in Fig. [2] as a function of Mu 3 ^. The notation for the 
three hierarchies is as in Fig. [TJ The thick lines are obtained using all available expansion 
terms whereas for the thin curves the highest order is set to zero. Thus the difference 
between the thick and the corresponding thin lines is a measure for the quality of the 
convergence. 

One observes a smiliar behaviour as at two-loop order: For small values of Mq 3> r, which 
correspond to small values of m^, both (h2) and (h3) provide good approximations. With 
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Figure 1: (a) mf , and 6 t obtained from the diagonalization of the top squark mass 

matrix as a function of the soft SUSY breaking parameter Mu 3) r. (b) (aj as a function 
of Mu 3 ,r using the parameters of Eq. (Ilip . The exact result is shown as solid black line. 

increasing Mu 3 ,r (h3) becomes worse whereas Qpp and Qpp for (h2) are still practically 
on top of each other. For values 300 GeV <^ M^ 3 ,r ^ 800 GeV the top squark masses are 
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Fi gure 2; (^ois as a function of Mh 3) r using the parameters of Eq. (jlip . Thick lines include 
all available terms whereas for the thin lines the highest terms are cancelled. 
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relatively close to each other which is the region of validity for (hi). For higher values 
one observes again a strong hierarchy between and m~ t2 and thus (h2) takes over. It 
is interesting to note that for each value of Mu 3t R there is at least one hierarchy with a 
small value of 5 app and thus an expected good approximation to the unknown exact result. 
Furthermore, the approximations show a significant overlap so that the whole range of 
Mu 3t n is covered. 

Let us in the following briefly discuss the numerical impact of the three-loop corrections 
computed in this paper. In Figs. [3] we show the strong coupling at the GUT scale, 
ai S( ^ CD ^(MGUT) with Mqut = 2 • 10 16 GeV as a function of the decoupling scale /id ec 
which is obtained by the following procedure. The starting point is ai^' MS (Mz)- In a 
first step we run in the SM from \i = M z to fi = yU dcc where the decoupling of the top 
quark and the SUSY particles is performed simultaneously and a^\ndcc) is transformed 
to ai SQCD ^(/idcc)- The use of the SQCD (3 function finally leads to «s CD ^(Mgut)- The 
thick lines in Fig. [3] correspond to this procedure, i.e., we use the following chain in order 
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Figure 3: ai S( ^ CD ^(MGUT) as a function of /idee- Thick and thin lines correspond to the one- 
and two-step scenario, respectively. Thin lines are only shown for three- and four-loop 
running. 

to arrive at Cd ^(Mgut) 

a f)m {Mz) run. a pVS M dec. a W*s» M run. a (SQCD) (MGuT) (13) 

For a degenerate supersymmetric mass spectrum the decoupling constant can be found 
in Eq. (E}. 

Alternatively, in order to obtain the thin lines we integrate out the top quark in a separate 
step with /i = M t (M t is the on-shell top quark mass) and transform afterwards ai 6 ^' MS 
to ai S< ^ CD ^(MGUT) hi analogy to Eq. (fTBl . Thus we have 



run. (5) i ms/ ]\/t \ dec. (6) 



MS/ 



onstant needed for the trans 
degenerate SUSY masses is given in Eq. ( JIUIl . 



GUT ) 



(14) 



The decoupling constant needed for the transition from ai 6 ^' MS to ai 3 ^ 00 "* in the limit of 
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In order to obtain the numerical results in Fig. [3] we have used the measured result for 
af\M z ) which reads [21] 

af(M z ) = 0.1184 ±0.0007, (15) 

Furthermore we have adopted a mSUGRA scenario with 

m = 700 GeV , m 1/2 = 600 GeV , tan/3 = 10 , A = , /i S usy > . (16) 

as input for sof tsusy [22] in order to compute the supersymmetric mass spectrum. Note 
that there is only a weak dependence of the general features of our numerical result on 
the particular spectrum. However, it is convenient to make use of a spectrum generator in 
order to obtain directly the DR values for the masses at the scale /idee- To our knowledge 
the running of the DR parameters is only implemented to two-loop accuracy which poses 
a slight inconsistency in our analysis. However, this is only an minor effect and does not 
influence the main conclusions. In order to get an impression about the numerical values 
for the physical masses we show the DR results for a typical scale /idee = 1000 GeV 

m t = 146.7 GeV , m h = 1022 GeV , m h = 1271 GeV , 

rriq = 1348 GeV , m~ g = 1326 GeV , 6 t = 1.26 . (17) 

At three-loop order the best appoximation is provided by the hierarchy (hi). In fact the 
quantity 5 app in Eq. (lT2~j) takes the value <5 app = 0.002. 

For consistency iV-loop running has to be accompanied with N — 1-loop decoupling rela- 
tions. Thus, we can show curves for N — 1,2,3 and 4 which corresponds to the (thick) 
dotted, dash-dotted, dashed and solid line, respectively. Within QCD the beta function 
is known to four-loop accuracy [23l[23], however, the supersymmetric analogue only to 
three loops jT6j[25l|26]£] As a consequence for the four-loop curve in Fig. [3] we only use 
three- loop running above /idee- 

/idee is an unphysical scale not predicted by theory. Thus, on general grounds, the depen- 
dence on /^dec has to diminish if higher order corrections are included. This is clearly visible 
in Fig. [3] where the dotted, dash-dotted, dashed and solid lines correspond to one-, two-, 
three- and four-loop running, respectively. Around the central scale of approximately 
1000 GeV all loop orders lead to predictions which are quite close. However, a variation 
of /idee leads to a relatively strong variation of the two-loop result which gets stabilized at 
three-loops and which furthermore gets to a large extend /idee independent at four loops. 
Actually, varying /j dcc between 100 GeV and 10 000 GeV changes ai S ^ CD ^(M GU T) by only 
0.07%. 

It is interesting to compare the variation of the individual curves with respect to /id ec 
with the experimental uncertainty induced from a.s' US (Mz) which is indicated by 
the band around the four-loop curve. The two-loop prediction is inside the band for 
300 GeV <^ /idee ^ 1800 GeV whereas the three-loop curve leaves the band only for 

1 The four-loop SQCD f3 function is not yet complete [27] . 
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/-tdec <: 13 000 GeV. It is also interesting to mention that all higher order corrections are 
very small for fi dec ps 650 GeV. 

Note that often /id ec = Mz is chosen for the matching between the SM and the MSSM. 
This choice leads to strong deviation at two-loops. At three-loop order the results are 
already quite stable which is further supported at four loops. 

Let us finally remark on the step-by-step decoupling of the top quark and the supersym- 
metric particles. The corresponding three- and four-loop results are shown as thin lines 
in Fig. El One observes even flatter curves than for the one-step scenario, however, the 
difference is numerically small and well within the uncertainty band. In this context we 
want to stress the wide range of //dec which is considered in Fig. |3j 



3 Low-energy theorem and Higgs-gluon coupling in 
supersymmetric QCD 

In Ref. [4] the following formula valid for all orders in perturbation theory has been derived 
in the framework of 

d = A? CD lnC« s (18) 

where 

describes the derivative with respect to the heavy mass m^. Thus the iV-loop corrections to 
C ag immediately leads to iV-loop corrections to C\. Since in Eq. (fl8"|) a logarithmic deriva- 
tive is taken and furthermore the dependence of ( as on rrih only occurs via ln(/^ 2 /m|) even 
the (N + 1) corrections of C\ can be computed once the renormalization scale dependence 
of ( as at (N + l)-loop order is re-constructed with the help of the renormalization group 
equations. 

The LET of Eq. ( TTB"]) can easily be extended to the case where more than one heavy quark 
is present. This version has been used in Ref. [7] in order to derive C\ for theories with 
several heavy quarks which couple in a Yukawa-like way to the Higgs boson. 

The extension of Eq. ( TTB"]) to NLO corrections in the framework of the MSSM has been 
considered in Ref. [TT]. Because of the different setup of our calculation, which is mainly 
due to the e scalars, we cannot take over the derivation of Ref. [TT]. However, following the 
same line of reasoning as in Ref. [I] we obtain a version of the LET which is appropriate 
for the decoupling constants computed in the previous chapter. For this purpose it is 
convenient to consider the bare decoupling constant (see Eq. (J4])) expressed in terms of 



3 Note the different normalization of the operator 0\ in Ref. 
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bare parameters. This leads to the LET in the form 

C? = D° h \n(° as . (20) 

D® contains derivatives with respect to bare parameters (indicated by the superscript 
"0") and can be written as 

A £ is the evanescent Higgs boson-e scalar coupling which is best defined through the 
corresponding part of the Lagrange density [T] 



1 

2 v " £j v° 



C £ = -^(M e °) 2 £ ^-^(A°) 2 s^. (22) 



For convenience we have also displayed the mass term for the e scalar. 
The derivative operators in Eq. ( I2~TT) are defined throughEl 
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where the prefactors in the top quark sector are obtained from the relations 

_ cos a 

sm p 

COS / \ 

V L * L = - 2m l^p + M z cos 2 ( 1 - g tan 2 W J sin(a + /3) , 
cos 4 

= -2m 2 —— + -M 2 Z sin 2 W sin(a + 0) , 
sm u o 

Kr = ^rl = 4^(-/i S usY sin a - A t cos a) , 
sm jD 
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with 



R(9 t 



cos 9 t — sin 
sin 6 t cos t 



All light quark masses are set to zero. Their averaged contribution is denoted by q and 
thus we have 

- 1_ ( ni + rn j n x -n t , 

m V 2 2 



3 In order to keep the notation simple we omit the superscript "0" in these expressions. 
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- _ 1_ / rii + nt a tm-th fi \ 
V 2 RR + 2 Ha j ' 

= y« =o 

K LR. R.L w ' 

= M 2 Z cos 2 W M - ~ tan 2 9 w ^j sin(a + /3) , 
V^ R = ^M|sin 2 ^sin(« + /3), 
V* L = M% cos 2 9 W ^-1 - ^ tan 2 6>y^ sin(a + /3) , 

= sin 2 W sin(a + /3) , 

where "it" and "cf ' denote generic up- and down-type squarks, respectively, and the labels 
tli = 5 and nt = 1 are kept arbitrary for convenience. with i, j = 1,2 are obtained in 
analogy to V^. 

After applying D® to of Section [2] we obtain the coefficient function C\ expressed 
in terms of bare parameters. Thus, in a next step one has to perform the parameter 
renormalizationQ Furthermore, it is necessary to take into account the operator renor- 
malization constant, often denoted by Z\\ [TS] . to obtain a finite result for the coefficient 
function which can then be compared to [T^[T5]. 

An alternative version of the LET (compared to Eq. ( l20l) ) is obtained by exploiting the 
fact that Z aa and Z a i g are independent of the parameters occurring in D®. Thus we can 
write 

d = D° h \n( as , (23) 

where it is still understood that C\ and ( as are expressed in terms of unrenormalized 
parameters. After computing C\ with the help of Eq. (123]) the parameters have to be 
renormalized as before, however, the operator renormalization constant is not necessary 
anymore. 

A third version of the LET reads 

Ci = DfclnCo.. (24) 

In this equation all quantities are expressed in terms of DR renormalized quantities and 
^sqcd^ exce pt the evanescent couplings (M e and A e ) which are renormalized to zero. It 
is very convenient to use Eq. (I24p since it directly leads to a finite result for C\. It is 
worth noting that the computation of C\ from Eq. (12^1) avoids the introduction of the 
evanescent coupling A e . This can be understood by considering the renormalized version 
of Dh in Eq. fl21|) where the last term vanishes due to the condition A 2 = (A°) 2 — 5A 2 = 0. 

4 The details are described in Ref. [HUTS] . 
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Due to the derivatives in Eq. (l2Tj) the expansion depth available for ( as is reduced. Nev- 
ertheless we can compare the results to the findings of Ref. [HJ[T5] where C\ has been 
computed from vertex diagrams. For all three hierarchies we found complete agreement 
for the first three terms in the mass difference, i.e. up to order {mf — m 2 ) 2 . Furthermore, 
for (hi) [(h3)] terms up to 1/?tisusy [1/^susy] com d De compared successfully and for 
(h2) all terms including 0(l/m~ ) and 0(l/m^ VSY ) agree. Thus the calculation of the 
decoupling constant together with the application of the LET provides an independent 
confirmation of the Higgs-gluon coupling at three-loop order. 

The LET in Eq. ( 120]) differs from the one presented in [11] by the term involving A° (see 
Eq. (I2ip ). Up to NLO it is possible to avoid such a contribution [TT][T3j. at three-loop 
order, however, a renormalization of the Higgs boson-e scalar coupling is mandatory (see 
Ref. [Hldn] for a detailed discussion) in case derivatives with respect to bare parameters 
are taken. 

4 Conclusions 

In this paper we have computed the three-loop SQCD corrections to the decoupling con- 
stant relating a s defined in full MSSM to the one defined in QCD. The occurring three-loop 
integrals have been evaluated by applying expansions in various hierarchies and thus re- 
sults are obtained which are valid in a large part of the parameter space. The decoupling 
constant constitutes an important ingredient in the relation of a s (Mz) and o^Mgut)- 
We have shown that the inclusion three-loop terms to the decoupling constant in combi- 
nation with four-loop corrections to the /3 function leads to results for « s (Mgut) which 
are practically independent of the decoupling scale /idee, where the effective theory is 
matched to the full one, even when considering a variation of /idee by more than two 
orders of magnitude. 

A further interesting application of the decoupling constant is its relation to the effec- 
tive Higgs-gluon coupling C\ which is obtained by simple derivatives with respect to the 
involved parameters. This calculation constitutes an independent check of the results ob- 
taines in Ref. [T4"j[T5] by an explicit calculation. In this paper we provide the corresponding 
LET which contains all features also present at higher orders in perturbation theory. It is 
valid to all orders in perturbation theory. We have checked that the renormalized version 
(cf. Eq. (|24p ) works including three-loop SQCD corrections. 
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Appendix: Exact one- and two-loop result for ( as 

In this Section we present the results for ( as up to two loops taking into account the exact 
dependence on the occurring masses. All parameters are renormalized in the DR scheme 
except M £ which is renormalized on-shell. 

In contrast to Eq. (jSJ) the coefficients of defined through 
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is presented. The results read 
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where Cp = 4/3, = 3, T F = 1/2, N t = 1, N q = 5, £ n is the Riemann zeta function, 
Z x = ln(/i 2 /m 2 ), L e = ln(/i 2 /M 2 ) and M e is the e scalar mass. Furthermore we have 
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where following abbreviations have been introduced 



A(x, y) =yj (1 - x - y) 2 - 4xy , 
Ch(x) =lm[Li 2 (e ix )\ , 

$i(x, y) =A~ 1 (x, y)|21n [|(1 + x-y- A(x, y))] In [1(1 -x + y- X(x, y))] + ±vr 2 
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The one-loop result agrees with Ref. [8] . The two-loop result has also been considered in Ref . [TO] , 
however, no compact result has been presented. Furthermore, the decoupling has only been 
considered within DRED, i.e., the transition from DREG to DRED has been performed in a 
separate step. 
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